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Abstract: We study non-perturbative aspects of the large ./V duality between Chern-Simons 
theory and topological strings, and we find a rich structure of large N phase transitions in the 
complex plane of the 't Hooft parameter. These transitions are due to large ./V instanton effects, 
and they can be regarded as a deformation of the Stokes phenomenon. Moreover, we show that, 
for generic values of the 't Hooft coupling, instanton effects are not exponentially suppressed 
at large N and they correct the genus expansion. This phenomenon was first discovered in the 
context of matrix models, and we interpret it as a generalization of the oscillatory asymptotics 
along anti-Stokes lines. In the string dual, the instanton effects can be interpreted as corrections 
to the saddle string geometry due to discretized neighbouring geometries. As a mathematical 
application, we obtain the 1/N asymptotics of the partition function of Chern-Simons theory 
on the lens space L(2, 1), and we test it numerically to high precision in order to exhibit the 
importance of instanton effects. 
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1. Introduction 

According to large iV dualities, the 1/N expansion of different observables in U (N) gauge theories 
can be reinterpreted as a genus expansion in an appropriate string theory. For example, the gauge 
theory free energy at finite volume has an expansion of the form 

oo 

F(g.,N)~'£g*>-*F g {t), (1.1) 

9=0 

where g s is the gauge theory coupling constant and t = g s N is the 't Hooft parameter. In large 
N dualities, g s is identified with the string coupling constant, the 't Hooft parameter becomes a 
geometric modulus of the string target space, and the amplitudes F g (t) are identified with free 
energies at genus g of a string theory. 

The 1/N expansion in gauge theory is an asymptotic expansion, and it is expected to have 
corrections of the form 

oo 

F(g s ,N) ~ Y,9s g - 2 F 9 (t) + 0(e- A(t ^). (1.2) 

9=0 
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These corrections are invisible in the traditional 1/N asymptotics 1 , and they are associated to 
large N instantons. Although large N instantons are typically built upon classical gauge theory 
instantons, they should not be confused with them. A classical instanton is a saddle point of 
the classical action, while a large N instanton is a saddle point of the large N effective action; 
see [1, 40] and [23] for explicit examples in the CP^ model and in two-dimensional Yang-Mills 
theory, respectively. In particular, the action A(t) of a large N instanton is a non-trivial function 
of the 't Hooft parameter which includes quantum corrections at all loops. In some cases, large 
N instantons can be interpreted in terms of D-branes of the string theory dual, as it happens for 
example in matrix models of non-critical strings. 

One could think that the instanton corrections in (1.2) are unimportant since they are 
exponentially suppressed at large N, but this is not always the case. As we move in parameter 
space they might become of order one and start contributing to the large N asymptotics. This 
scenario was advocated long ago by Neuberger [41] in order to explain the occurrence of critical 
points in Hermitian matrix models and the third-order large N phase transitions of unitary 
matrix models [24, 47]. The same mechanism explains as well the Douglas-Kazakov transition in 
two-dimensional Yang-Mills on the sphere [16, 23]. In these instanton-driven phase transitions, 
as we change the 't Hooft parameter at fixed g s , the action of the instanton A(t) vanish at some 
finite value t = t c . This leads in general to different 1/N expansions for the regions t > t c and 
t < t c . 

On the other hand, if we regard g s as a complex variable, we should expect that the 1/N 
asymptotics changes discontinuously as we change the argument of g s (and keep |i| fixed). In 
classical asymptotic analysis, this discontinuous change is the well-known Stokes phenomenon. 
The reason for this phenomenon is in fact the same one that underlies instanton-induced phase 
transitions: an exponentially small quantity, due to a sub-leading saddle point, becomes less 
and less suppressed as we change the argument of the expansion variable. Along the so-called 
anti-Stokes line, the contribution of this saddle is of order one and has to be included in the 
asymptotics. 

In general, a large N gauge theory with instanton sectors, should display a complex pattern 
of large N phase transitions, combining the instanton-driven phase transitions at fixed g s with 
the Stokes- like transitions as we vary the argument of g s . If the gauge theory that we are studying 
has a string dual in at least one of the phases, it is certainly interesting to understand what are 
the implications of these transitions in the string picture. 

In this paper we will analyze these issues in a family of models which have well-understood 
large N duals, namely Chern-Simons theory on the lens space L(p,l). This gauge theory is 
described by topological string theory on a Calabi-Yau space given by an A p fibration over a 
two-sphere [22, 3, 27, 26]. However, this large N string dual is intrinsically semiclassical, since 
it describes a generic but fixed saddle-point of the gauge theory. A full non-perturbative study 
of this duality, where we sum over all saddle-points, reveals a surprisingly rich phase structure 
as we move in the complex plane of the 't Hooft coupling. The main results of our analysis are 
the following: 

1. For each complex value of t, the large N asymptotics of the non-perturbative free energy 
is dominated by a fixed string target geometry, i.e. by a point in the Calabi-Yau moduli 
space. However, this "saddle geometry" might change as we vary t. 

1 By traditional or classical asymptotics, we mean asymptotic expansions in which non-analytic terms are not 
taken into account. There are refinements of traditional asymptotics, which go sometimes under the name of 
"hyperasymptotics," where these terms are included in a systematic way, see [7] for an overview. 
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2. For some values of t, the 1/N asymptotics is given by a conventional genus expansion. 
There are non-perturbative effects due to the contributions of neighbouring geometries, 
but they are exponentially suppressed. For other values of t, the asymptotics has an 
oscillatory behavior due to large N instantons. This type of behavior was first observed 
in the context of matrix models in [10], and further studied in [17, 18]. In particular, the 
large N asymptotics is no longer given by a genus expansion around a fixed geometry, and 
corrections due to neighboring geometries are crucial, already at the next-to-leading order. 

3. The change of saddle geometry as we change the complex 't Hooft parameter can be re- 
garded as a generalized or "deformed" Stokes phenomenon. In the limit of vanishing 't 
Hooft coupling, the asymptotics changes discontinuously as we change the argument of the 
string coupling constant, and we recover the classical Stokes phenomenon. For finite t the 
phenomenon is smoothed out, and the asymptotics changes continuously. 

In practice we have focused on the simplest, nontrivial model in the family, namely L(2, 1) = 
MP 3 . The gauge theory saddle points are characterized by a symmetry breaking pattern U (N) — > 
U(Ni) x U(N2). The dual Calabi-Yau manifold is the so-called local Fo geometry, where Fo = 
P 1 x P 1 is a Hirzebruch surface. It has two Kahler parameters, corresponding at large radius to 
the sizes of the two P 1 s in Fo- The partial 't Hooft couplings t\, <2, where U = g s Ni, parametrize 
the Kahler moduli space near the orbifold point described in [3], and their sum t\ + ti is fixed 
to be the total 't Hooft parameter of the gauge theory. Each saddle geometry is characterized 
by an "equilibrium value" for ti- The phase diagram in the region \t\ > 0, < Arg(<7 s ) < tt/2 is 
shown in Fig. 7. 

Our analysis shows very clearly that the large N instantons analyzed in [10, 35, 36, 34, 17, 
18, 43] are crucial in order to understand large N dualities. To be precise, for generic values 
of the 't Hooft parameter, the correct large N asymptotics of the exact gauge theory partition 
function involves these instanton effects already at next-to-leading order. In order to dispel any 
doubt about this, we study in detail the partition function of Chern-Simons theory on MP 3 , 
for imaginary values of the coupling constant. These are precisely the the "physical" or "on- 
shell" values of the coupling constant in Chern-Simons theory. Our analysis shows that these 
values correspond to an anti-Stokes line with oscillatory asymptotics, in which the traditional 
genus expansion has to be supplemented with instanton corrections. We test numerically our 
predictions for the instanton corrected asymptotics against the exact non-perturbative Chern- 
Simons partition function, and we obtain an impressive agreement. 

Our analysis also clarifies the issue of background dependence in topological string theory. 
Let us focus again on the example of Chern-Simons theory on MP 3 . Perturbative topological 
string theory depends on a choice of background (t\,ti). The gauge theory partition function 
is background independent, in the sense that it only depends on the total 't Hooft parameter 
t = t±+t2 (which is, as needed for consistency, a modular invariant in the dual topological string). 
However, the large N asymptotics of the partition function is "peaked" around a particular 
value of (ti, £2) which depends on t and corresponds to a particular saddle geometry. This fixed 
background emerges as an equilibrium value as a consequence of the large N limit, but it has no 
meaning at finite N. In fact, in the full theory, the value of the background (specified by, say, 
the value of ti) is not fixed, since ti is rather an internal or fluctuating variable which we have 
to sum over. But as usual in statistical mechanics, in the thermodynamic limit N — > 00 this 
fluctuating variable is peaked around an equilibrium value. 

We would like to point out that the role of instanton effects in the 1 /N expansion and their 
string theory interpretation has been addressed before in different contexts. For example, the 



- 3 - 



instanton-induced Douglas-Kazakov phase transition has been interpreted as a breakdown of 
the large N string dual [16]. In AdS/CFT at finite temperature, an incarnation of the Gross- 
Witten-Wadia phase transition has been argued to correspond to a breakdown of the geometric 
description of the string target [4] 2 . The importance of exponentially small effects, and their 
effects on the string target space, has been emphasized in [37], which studies as well the Stokes 
phenomenon in the open string moduli space. In [42], a large N phase transition in N = 4 super 
Yang-Mills theory on K3 was found, separating a phase dominated by the zero-instanton sector 
from a phase dominated by an instanton sector with large instanton number. In the context 
of matrix models, general techniques to analyze large N phase transitions have been developed 
along the years, culminating in the formulation in terms of Boutroux spectral curves in [9]. Of 
course, the fact that the genus expansion can not capture the large N asymptotics was discovered 
in the context of matrix models and explained in detail in [10]. 

This paper is organized as follows: In Section 2 we review the structure of instantons correc- 
tions in matrix models and topological strings. In Section 3 we present our model, the Chern- 
Simons theory on the lens space L(p, 1), and its matrix model realization. We then specialize 
to the p = 2 case and review the large N duality between the perturbative expansion of the 
matrix model in a fixed filling fraction (Ni,N2) configuration and the topological string on the 
local P 1 x P 1 geometry, in a fixed background t\ , . In Section 4 we study the large N phase 
diagram. We start by considering two particular cases: real and imaginary g s , which correspond 
to two phases of the gauge theory with two distinct dominant saddle configurations, each dual 
to a fixed background geometry. By following the analysis of Section 2, we explain the structure 
of the large N asymptotic of the two phases and provide very precise numerical checks of our 
predictions. We then complete the study of the phase diagram for generic complex g s and inter- 
pret the large N phase transitions between the various phases as due to the Stokes phenomen. 
We end this section by commenting on the issue of the background independence. In Section 5 
we discuss the phase diagram of the cubic matrix model. Thanks to its close relation with the 
Airy function, the prototypical example of the Stokes phenomenon, this example clarifies the 
nature of the large N phase transitions between dominant saddles/backgrounds which we found 
in Chern-Simons theory. Section 6 contains conclusions and a list of open issues. Finally in the 
appendix we collect some details of the model for different submanifolds of moduli space. 

2. Instanton corrections in matrix models and topological strings 

In this section we review instanton corrections in matrix models/topological strings, following 
mainly [10, 36, 18]. 

Multi-cut matrix models and topological strings are characterized, at the perturbative level, 
by genus g amplitudes of the form 

F g (t) = F g (*!,-.. ,t s ), (2.1) 
where, in the case of matrix models, U are partial 't Hooft couplings 

ti = g a Ni, i = l, (2.2) 

2 The Hawking-Page transition in AdS/CFT is a large N phase transition leading to a change of topology, but 
it involves effects of order exp(— TV 2 ) and it is not driven by instantons. 
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Ni is the number of eigenvalues in the i-th cut and s is the number of cuts. The perturbative 
partition function is given, as usual, by 

Z V {N U --- ,iV s ) = exp[f> a (t)^- 2 ]. (2.3) 

A choice of moduli or 't Hooft parameters t = (ti, ■ ■ ■ , t n ) in this partition function is called a 
choice of background. In matrix models and in gauge theory, the quantity 

n 

t = J2 t i = 9sN (2.4) 

i=l 

where TV" is the rank of the U(N) gauge group, is fixed. Background-independent quantities 
should only depend on g s and N, while the perturbative partition function Z p depends on a 
choice of t{ and is therefore a background dependent quantity. 

The non-perturbative partition function Z^ p was introduced in the context of matrix models 
with a polynomial potential V(z) in [13, 10, 17], and it depends on a choice of contour C in the 
complex plane. This contour is the integration path for the eigenvalues of the matrix integral. 
We then have 

zU^9s) = n^ ( A)e-^^ (2.5) 

' Jc i=i 

where A(A) is the Vandermonde determinant. Of course, the contour C is chosen in such a way 
that this integral converges. By standard saddle-point techniques, we can always deform the 
path C into a sum of paths Cj, which go through the critical points of V(z) and are paths of 
steepest descent [20], 

s 

C = J2(kC k , (2.6) 
fc=i 

where we have assumed that V(z) has s critical points. Therefore 

Z%(N,g s )= Y, C? 1 ~-C?'W, -,*.)■ (2-7) 

Nl +-+N S =N 

where 

^- k«, ' ' • /<•«. S £ A2(A) w ™ 

To write this equation, we have split the iV eigenvalues into s sets of Nj eigenvalues, I = 1, . . . ,s, 
which are denoted by 

{^}k I= l,...,N n 1 = 1,..., 8. (2.9) 

Each of the integrals (2.8) has a small g s asymptotic expansion given by (2.3). 

Equation (2.7) expresses a background-independent quantity Z% p as a sum of background- 
dependent quantities Z p (t). This sum can be in turn regarded as a sum over matrix model 
instantons, which have been identified long ago in terms of eigenvalue tunneling [13, 46] . Formally, 
we can write the sum (2.7) as the perturbative partition function coming from a fixed background 
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Ni = N* , plus an infinite sums of corrections for the remaining values of TVj. These corrections 
are non-perturbative in g s . We write (2.7), schematically, as 

Z%{N,g s ) = Z p (N* ir -- ,N;) + 0(e~ 1 ^) (2.10) 

In order to be more concrete, we will restrict ourselves to models where s = 2, or equivalently, 
models with two moduli t\, t2- Up to an overall normalization, we can set (l = 1, £2 = £. 
The detailed form of (2.10) depends on the choice of background. If the background is on the 
boundary, we have (N^N^) = (N, 0) or (0,7V). Otherwise, we say that the background is an 
interior point. The expression for (2.10) when the background is on the boundary was worked 
out in [36]. Assuming for concreteness that the boundary is at (N, 0), we have 

l 2 12 



1 + E(|^72 G ^ + 1K^ exp 



£A\ 
9s J 



E i (2g i +m i -2) 

V V V ^ F^ p^t-fS^ih 

k\m 1 \...m k \ 91 "' 9k 1 ; 



k mi >o Si> i_!^ 



(2.11) 



In this equation we have introduced the following notations. Z p (N) is the partition function for 
the one-cut model where all eigenvalues sit in the first critical point. The functions F g (ti,t2) are 
defined by 

Fgituh) = Ff(t 2 ) + F g {h,h), (2-12) 
where F g (t) are the genus g free energies of the gauged Gaussian matrix model with 't Hooft 
parameter t. Fg m ^ denotes the m-th derivative of F g w.r.t. s, which is defined as 

s = \(h-t2). (2.13) 

All derivatives are evaluated at t\ = t and t 2 = 0. G 2 (£ + 1) is the Barnes function 

e-i 

G 2 (£+l) = [J n! ( 2 - 14 ) 

n=0 

Finally, 

A(t) = d s F and g = exp(d 2 F ). (2.15) 
At leading order in g s , we have 

V ™- ... m fc ! 31 " " 9k 1 ; 

k ™.>o 3l >i-f (2.16) 

£ 3 



= 1 - g s (£d s F 1 (t) + - d 3 s F (t)) + 0(<? s 2 ). 



The expansion around an interior point was studied in [10, 17, 18]. It is given by 



k mi >o 9i>1 _^ fc!mi! '--^ ! 91 

= z p (NZ,NZ)le^ + gs (e'^F{ + ^e"„K') + o( g 2 s )\. 



(2.17) 
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The derivatives of the free energies F g are again w.r.t. s defined in (2.13). The theta function 
®u,u with characteristics (fi, v) is defined by 

@Hv{u, T) = ^e*"^ - ^" e wi(n+ii-Ne)T(n+n-Ne) e 2m(n+n)v (2.18) 

and it is evaluated at 

u = r = —Fq. (2.19) 
In the above equation, we have denoted 

N* 

e=-^, C = e 2i ^ (2.20) 

in order to make contact with the notations of [18]. In (2.17), the derivatives of the theta function 
(2.18) are w.r.t. u, therefore each derivative introduces a factor of n + \x — Ne in the sum (2.18). 
Notice that, as emphasized in [10], the expression (2.17) gives a large N asymptotics which goes 
beyond the genus expansion. In particular, the leading terms in the free energy are of the form, 

F = g; 2 F (N*) + F^N*) + log 9^ + • • • , (2.21) 

and the theta function leads typically to an oscillatory large N asymptotics. 

The above expressions, (2.11) and (2.17), are for the moment being merely formal. In 
particular, it is not clear if they provide reasonable asymptotic expansions of the original, non- 
perturbative partition function. There are again two different cases, applying to (2.11) and 
(2.17): 

1. The expression (2.11) gives an admissible asymptotic expansion if 

Ref^)>0. (2.22) 

In this case, the classical, small g s asymptotic expansion is given by the genus expansion of 
Z p (N). The terms with i > 1 give exponentially suppressed corrections to the asymptotics. 

2. The expansion around a generic point N* in (2.17) is admissible if the filling fractions, as 
determined from the spectral curve, are real, i.e. 

N* 

e t = -^r G R, (2.23) 

and if 

Re f F MH) = 0. (2.24) 



s 



The conditions (2.23), (2.24) were first spelled out in detail in [13, 14]. As noticed in [14], 
they can be written as 

lJ-Ip T y(x))=0, (2.25) 



t 7 7 2vri' 



where y(x) is the spectral curve of the matrix model, and 7 is any cycle on it. In writing these 
equations, we have followed the conventions of [32], so that 

y\x) = V'(x) 2 + ■■■ . (2.26) 
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A curve with the property (2.25) is called a Boutroux curve. The condition (2.24) makes sure 
that the term involving u in the theta function (2.18) is oscillatory. This in turn guarantees that 
the sum over n will be peaked around 

n + ^-iVewO (2.27) 

so that the derivatives of the theta function involved in (2.17) are of order one, see [19] for a 
related discussion. The condition that the curve is Boutroux gives a nice geometrization of the 
search for good asymptotic expansions in matrix model theory, and it has been recently developed 
in much detail in [9]. 

The above conditions imply that, in order to have a well-defined large ./V asymptotics, the 
backgroud can not be arbitrary. As we change the parameters of the theory, this "dominant 
background" will change. In particular, it can change from a boundary point to an interior 
point. We will see examples of this in the Chern-Simons model discussed in this paper. 

The search for a dominant background is however subtle, since as it is well-known, the 
perturbative amplitudes F g (ti) have a finite radius of convergence. The full space of 't Hooft 
parameters has a nontrivial global structure. As we reach the boundary of the convegence region, 
we have to perform a duality transformation to a different frame. In particular, the condition 
(2.24) is attached to a particular frame and region, in which Fo is well-defined. One could worry 
that the search for dominant backgrounds is not a well-posed problem, globally. However, we 
will now argue that this is not the case for two reasons. The first one is that the condition on the 
curve being Boutroux is invariant under duality transformations. Second, as shown in [18], each 
term in the l/N expansion of (2.17) is invariant under duality transformations, up to a phase. 
This means, in particular, that 

(2.28) 



p g7 2 F +F lf) 
e ^M,^ 



is invariant and it is therefore a well-defined quantity in the global moduli space. Looking for 
the maxima of this function on the subspace of Boutroux curves in moduli space is therefore 
a well-posed global problem, i.e. independent of the duality frame. In each region of moduli 
space in which Fo, Fx, Mj!/ are not singular, the maxima at large N satisfy the condition (2.24). 
Therefore, the problem of solving (2.23), (2.24) can be lifted globally in a consistent way. 

The structure of instanton corrections that we have reviewed here was originally derived in 
the context of matrix models. However, it was proposed in [35, 34, 18] that, since it only de- 
pends formally on the data of the spectral curve, these corrections should be present in topological 
strings on local geometries, including toric geometries. In this sense, a suitable example is pro- 
vided by topological string theory on A, p fibrations on P 1 . These models have a non-perturbative 
definition in terms of U(N) Chern-Simons gauge theory on a lens space [3], and they can be 
explicitly reformulated as a multi-cut matrix integral. They constitute a privileged arena for 
constructing the large N asymptotics presented above and for studying the role of instanton 
corrections. In the next section we present what is known about these models from the point of 
view of Chern-Simons theory, matrix models, and topological string theory. 



3. Large N duality for Chern-Simons theory 

In this paper we will study Chern-Simons theory on the lens space L(p, 1). This model has an 
open string description as an A-type topological string defined on a Calabi-Yau threefold given 
by the cotangent bundle over the lens space T*L(p, 1). This is just a particular example of the 
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more general realization, due to Witten [49], of Chern-Simons theory on a three-manifold M as 
a type A topological string on T*M. 

For M = S 3 , which corresponds to the lens space with p = 1, the Gopakumar and Vafa large 
N duality [22] provides a description of U (N) Chern-Simons gauge theory in terms of closed 
topological string theory on the resolved conifold 0{— 1) + 0{— 1) — > P 1 . This duality can be 
generalized [3] by performing a Z p quotient on both sides of the Gopakumar- Vafa duality. This 
leads to a duality between Chern-Simons theory on S^/Zp = L(p,l) and the closed topological 
A-model string on an A p fibration over P 1 . This target space, unlike the resolved conifold, has a 
non trivial Kahler moduli space and besides the large radius phase it has non-geometric phases. 
In particular, the large N duality relates the perturbative regime of the Cherns-Simons gauge 
theory — small 't Hooft couplings — to the orbifold regime of the string theory. This duality has 
been tested in [3] by comparing the perturbative expansion of Chern-Simons theory around a 
fixed flat connection, to the orbifold topological strings amplitudes computed by mirror sym- 
metry. Further tests of the duality have been done in [27, 26]. In the following we will review 
some of these results. We will however go beyond the perturbative side of the duality and we 
will study the exact, non-perturbative gauge theory partition function. As we will see, in the full 
theory, the choice of filling fractions is not arbitrary, but it is dictated by the only parameters of 
the theory, the Chern-Simons level k and the rank of the gauge group N. 

3.1 Chern-Simons theory on lens spaces 

The lens space L(p,l) is a three-manifold that can be obtained by gluing two solid 2-tori along 
their boundaries after performing the SL(2,Z) transformation, 

(,„ 

This description makes it possible to calculate the partition function of Chern-Simons theory 
on these spaces [48]. SL(2,Z) transformations lift to operators acting on 7i(T 2 ), the Hilbert 
space obtained by canonical quantization of the Chern-Simons theory on the two-torus. This 
space is the space of integrable representations of a Wess-Zumino- Witten (WZW) model with 
gauge group G at level k, where G and k are respectively the Chern-Simons gauge group and 
the quantized coupling constant. In particular the partition function of the U(N) Chern-Simons 
theory on the lens space L(p, 1) is given by: 

Z(L(p,l)) = (0\U p \0), (3.2) 

where U p is the lift of (3.1) to an operator on li(Y 2 ) and |0) is the vacuum state, corresponding 
to the trivial representation. This partition function has the following explicit form (up to an 
overall constant) [44]: 



2nl 



(_1)|A+I e k p C fi7r 

neZ N /pZ N wew 



^e(w)exp{^(p 2 -2p-(kn + w(p)) + (kn + w(p)) 2 )y (3.3) 



( \pk) N / 2 n(z7 N /r , T N wG W kp 



In this expression, k = k + y, where y denotes the dual Coxter number, |A+| is the number of 
positive roots of G, and p is the Weyl vector. Using Weyl's formula 

e{w)e w(p) = JJ 2 sinh f ( 3 - 4 ) 
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we can also write the partition function as: 



i|A+| 



Z^-g— E ^~" P - n n 2sin(^(Mn,-n,) + i-,)). (3.5) 

I '^"V n ^7N /„jn l<i<j<N W 



This expression for the partition function can be understood as a sum over saddle-points. For 
U (N) Chern-Simons theory on L(p, 1) = §3/Z p , these are flat connections, which can be obtained 
by considering embeddings of the first fundamental group into U(N), modulo gauge transfor- 
mations. Since iri (L(p, 1)) = Z p these embeddings are given by iV-component vectors n, whose 
entries take values in Z p . However, since the residual Weyl symmetry Sn of the U (N) gauge 
group permutes the different components of n, the flat connections are rather labeled by parti- 
tions Nj, where Nj is the number of entries in n which is equal to j — 1. The Nj label the choice 
of vacuum 

U(N) -»• U(Ni) x • • • x U(N P ). (3.6) 
It is then possible [3] to rewrite the partition function as a sum of contributions of flat connections: 

z = E z w ( 3 - 7 ) 

{Nj} 

where 

Z (^)= /j 1 D/' 2 n V» £ e(w)e W {^(p 2 -2p.(kw'(n) + w(p)) + (kw'(n) + w(p)) 2 )}. 
{ lple) llj j- WjW , eW k,p 

(3.8) 

In this expression, the vector n is any vector with Nj entries equal to j — 1. The sum over Weyl 
permutations w 1 6 W guarantees that the resulting object is gauge- invariant. As shown in [31] 
it is possible to rewrite (3.8) matrix integral 



W = " n / n§e-4&— ™» «=n(2smh^) 2 , (3.9) 

^1 ^ 1=1 l<j 



where 



9s = 2 4 (3-10) 



and the pre-factor involves the gauge-theory action of each flat connection: 

9s V"9s 



It can be easily checked that, although the expression (3.9) involves a choice of vector n, any two 
choices related by a Weyl permutation lead to the same matrix integral. The expression (3.9) 
can be regarded as ap-cut matrix model, where N^ eigenvalues sit around the point 2iri(k — l)/p, 
i.e. it is a matrix model with fixed filling fractions. It has in particular an asymptotic large N 
expansion of the form 

oo 

F Nk = log Z Nk = E F g {ti)9 2 r 2 (3-12) 

9=0 
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where 

U = g s N t (3.13) 

are the partial 't Hooft parameters. In the case p = 2, calculations around the Gaussian point 
of the matrix model give the genus zero free energy 

F (t 1 ,t 2 ) = -—t 2 + log(-4)M 2 + F G (ti) + F G (t 2 ) + F p (t!, t 2 ). (3.14) 

In this equation, the first term comes from the instanton action, the second term comes from 
the overall measure of the matrix integral, and F^(t) is the Gaussian matrix model genus zero 
amplitude, 

Ff(t)= 1 -t*(logt-l). (3.15) 
Finally, ^(ii,^) comes from fatgraphs of genus zero, and it is given by 



F$(h,t 2 ) = ^(4 + 6t?t 2 + m\t 2 + 6ht 3 2 + 4) 

' - (4tf + 45tf t 2 + 22ht\t\ + 1500^ + Tlht\t\ + 45M| + 4t|) + 



(3.16) 



345600 

For the genus one term one has: 

i ? i(ti,t 2 ) = if(iVi) + if(iV 2 ) + FP(ti,* 2 ) (3.17) 

where 

F 1 G (iV) = C / (-l) + ^log(iV) (3.18) 
is the Gaussian matrix model contribution, and 

Ff(ti,t 2 ) = -^g(t? - 6ti< 2 + i|) + ■ ■ ■ (3.19) 
comes from fatgraphs of genus one. Higher genus free energies can be computed analogously. 
3.2 The dual topological string 

We will now restrict ourselves to p = 2. The large iV dual model of Chern-Simons theory on 
L(2, 1) is the topological string on the anti-canonical bundle of the Hirzebruch surface Fq. This 
geometry has two P^s in H 2 (Y) and a compact four cycle in H^Y). We denote their associated 
classes by A±,A 2 and B respectively. The class C = A 2 — A\ does not have a dual cycle in 
H^iY) and thus it will correspond to a non-normalizable modulus of the theory; one can then 
take A = A 2 as the normalizable modulus. 

The mirror geometry is encoded in a family of elliptic curves S, which can be written as 

Z\X 2 + X + 1 - + X + Z\X 2 ) 2 - Az 2 X 2 

V — 2 ' (3.20) 

The periods of the one form A = log ydx / x are given by the solutions of the associated Picard- 
Fuchs system: 

d = 22(1 - 4z 2 )g - Aziil - 8z lZ2 ^ 2 - 6zi& + (1 - 6z 2 )6, 

(3.21) 

C 2 = zi(l - 4*i)f? - AzU 2 2 - 8z lZ2 ^ 2 - 6z 2 & + (1 - 6*i)£i, 
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where 

«■ " k (3 ' 22) 

At the large radius point z\ = z 2 = 0, the Picard-Fuchs system has two single logarithmic 
solutions, which give the mirror maps at large radius 

Ti = - log z x - 2zi - 2z 2 - 3zf - \2z x z 2 - 3z\ H , 

(3.23) 

T 2 = - log z 2 - 2zi - 2z 2 - 3zj - \2z x z 2 — 3z\ H . 

As expected, the quantity 

T 1 -T 2 = log-, (3.24) 

^i 

corresponding to the class C, does not receive any instanton corrections. It only depends on the 
"bare parameters" of the model z\ , z 2 and has to be regarded as a ■parameter. As we will see in 
a moment, it has a very natural meaning in the dual Chern-Simons theory. 

In the moduli space of this model there is another point, discovered in [3], which makes 
contact with Chern-Simons perturbation theory. This point is called the orbifold point, and it is 
defined as the point x\ = x 2 = in terms of the variables: 

si = 1 - — , x 2 = 7 1 v . (3.25) 



Z 2 . fZ2 



0-s) 



The periods near the orbifold point can be obtained again as solutions to the Picard-Fuchs 
system. They have the structure, 

(71 = - log(l - Xl) = ^ CmfiX'T, 
m 

C"2 = y ] C m> nXi X 2 , ^ 
m,n 

F a2 = cr 2 log(a;i) + ^dm^x^x^ , 

m,n 

where the coefficients c m ,n and d m ,n are determined by the following recursions relations [3] 

(n + 2 - 2m) 2 

C m n —Cm—l,n~T~? \7 TT' 

4(m — n)(m — 1) 

1 



n(n — 1) 



(cm,n-2(n -m - l)(n - m - 2) - c m _i in _2(n - m - 1) ), 



_4-i,n(n + 2 - 2m) 2 + 4(n + 1 - 2m)c m) „ + 4(2m - n - 2) (3.27) 

d m ,n — . i \ , -. n > 

4(m — n)(m — 1) 

4,n = 7 1 ^ {d m ,n-2{n - m - l)(n - m - 2) - d m _i n _ 2 (n - m - l) 2 
n(n — 1) 

+ (2n - 2 - 2m)c rn _i i „_ 2 + (2m + 3 - 2n)c mi „_ 2 ). 

The orbifold point gives the large iV solution of Chern-Simons theory on L(2, 1), around a fixed, 
arbitray flat connection, in the sense that the perturbative topological string amplitudes in this 
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frame, Fg (<r i , 02) > correspond to the amplitudes in the 1/N expansion (3.12). The dictionary 
relating the variables is 

h = -(cri+a 2 ), t 2 =-(0-1-0-2)- (3-28) 
Notice that the parameter of the theory 

log-=<7i = 2(ti+t 2 ), (3.29) 

is proportional to the total 't Hooft parameter of Chern-Simons theory g s N. This is in fact 
required by the large N duality: although t\ 2, the partial 't Hooft parameters, change in a 
natural way under the action of the symplectic group, the total 't Hooft parameter t must be 
a symplectic invariant, since it is a parameter of the gauge theory. The derivative of the genus 
zero free energy of Chern-Simons theory w.r.t. the variable s defined in (2.13) is given by the 
following combination of the periods: 

^ = y + ^ 2 -log(-4)0 2 . (3.30) 

By integrating w.r.t. a 2 and expressing the result in terms of t±,t 2 one recovers the genus zero 
free energy of the matrix model expansion. 

In the following it will be useful to express some quantities in terms of modular forms. It 
turns out that [2, 25] 

1 d 2 F 

T = w (3 - 31 » 

is indeed a modular parameter for the curve (3.20), and in particular Imr > 0. The modulus 

1 oo.l/ « 1 



u 



x\xl + - (VT^~ 1 + n , (3.32) 



8^1^ 1 ^ 2 V v 7 !^!, 
is closely related to the modulus of a Seiberg-Witten curve [45] with 

T SW = 2r + 1, (3.33) 
and it can be expressed in terms of elliptic theta functions: 

1? 4 - -i? 4 

u = 4 2 (2t) = 1 - 32g + 256g 2 + • • • (3.34) 

^3 

where q = e 2ntT . Notice that the orbifold point in the moduli space of local Fo corresponds to 
u = 1, i.e. to the monopole point of the Seiberg-Witten curve. 

Finally one can check that the genus one free energy is given, in terms of modular forms, by 

Fi = - log iKrgw) + 2C'(-1) + I log I , (3.35) 

In particular, we find 

Ff(h,t 2 ) = -\ogr)(T SW ) + ^log . (3.36) 

The partition function of Chern-Simons theory on L(2, 1) can be written in terms of a two- 
cut matrix model [31, 3]. The cuts are centered around x = 0, x = ni, and they are precisely the 
cuts which appear in the mirror curve (3.20). If we write them as 

(—a, a), (7ri — b, iri + b), (3.37) 
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in terms of the endpoints a, b, it is easy to see that they are determined by the equations (see 
also [27]) 

cosh a = —— + ,/— = l + 2ti + ti(ti+* 2 ) + -- - , 

2-v/iT V z i , 

1 /F" (3 ' 38) 

cosh b = - — — + . — = 1 + 2t 2 + t 2 (*i + h) + ■ ■ ■ ■ 

3.3 Two special slices 

The above results simplify considerably along particular submanifolds in moduli space, which 
will be relevant in the analysis of the phase structure in the next section. 

The slice t 2 = corresponds, from the point of view of the gauge theory, to an expansion 
around the trivial flat connection. Up to a trivial rescaling of the coupling constant g s , the theory 
reduces to Chern-Simons theory on the three-sphere S 3 . From the point of view of topological 
string theory, the theory is equivalent to the resolved conifold. The genus zero and genus one 
free energies are given by, 

if W = ^-^-Li 3 (e-<) + C(3), 

U b , (3.39) 

Ff(t) = -± - 1 log(l - e-«) + 1 log* + C'(-l) - y 2 lo % N - 

The instanton action (2.15) on t\ = t, t 2 = can be evaluated in closed form (see the Appendix) 
and it is given by 

A{t) = 2 Li 2 (e" i/2 ) - 2 Li 2 (-e-* /2 ). (3.40) 

The other relevant slice, for us, corresponds to the symmetric slice t\ = t 2 . In terms of the 
orbifold coordinates, this slice is parametrized by 

X2 = a 2 = 0, X! = l- e- ai . (3.41) 

On this slice we have the following exact expressions for the prepotential and for the modular 
parameter in terms of the flat coordinate o\ : 

where K{k) is the standard elliptic integral. The proof of the above identities is given in the 
Appendix. 



4. Large N phase diagram 

So far we have just recalled the perturbative large N duality between Chern-Simons theory on 
L(2, 1) expanded around a particular flat connection (corresponding to a choice of filling fractions 
Ni,N 2 ), and topological string theory on local Fo with Kahler parameters t±,t 2 . 

However, as we explained in the introduction, when studying the large N expansion of 
the non-perturbative partition function, the choice of filling fractions (dominant saddle) is not 
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arbitrary but it is dictated by the value of the total 't Hooft coupling t. We will now study the 
large N phase diagram for the Chern-Simons theory on L(2, 1), i.e. we will determine the saddle 
which gives the dominant contribution to the non-perturbative partition function as we move on 
the complex plane of t = g s N (equivalently, on the complex plane of g s ). 

We first notice two exact symmetries of the non-perturbative partition function. If we 
conjugate g s , we find 

Z(N,g* s ) = (-iy^Z*(N,g s ), (4.1) 

where c(N) is a half-integer that only depends on N. This symmetry is valid for all p. The 
second symmetry is related to changing g s — > —g s . If we write 

N 

Z(N,g s )=Y J Z(N,N 2 ,g s ) (4.2) 

we have that 

Z(N,N 2 ,-g s ) = (-ly'^e-^i-l^-^ZiN^ - N 2 ,g s ), (4.3) 

where c'(N) is again a half- integer that only depends on N. 

In determining the phase structure, we are interested in knowing which is the N 2 which 
contributes the most to the non-perturbative partition function (in absolute value). As we will 
see, for N large, the sum in (4.2) will be very peaked around a particular saddle value (N^N^). 
We see, from (4.1), that conjugating g s does not change this saddle, while from (4.3) we see that 
the phase diagram is invariant under the simultaneous change 

g s ^-gs, (NZ,NS)+>(NZ,NZ). (4.4) 

In practice, this means that in order to determine the phase diagram, we can restrict ourselves 
to the first quadrant of the complex g s plane 

< Arg( 5s ) < |. (4.5) 

In order to proceed, we first study the limiting cases of real, positive g s and purely imaginary g s . 

4.1 Real g s : the boundary expansion 

The classical action of the flat connection (3.11) has the value 

S ills t(N 2 ) = yiV 2 . (4.6) 

When g s is real and positive, it gives an exponential suppression for all the non-trivial flat 
connections N 2 > 0. Therefore, at least when the 't Hooft coupling is sufficiently small, the 
dominant saddle point is expected to be at the boundary (N^N^) = (N, 0). In fact, it is easy 
to see that for any real, positive t, this boundary point is the dominant saddle. To prove it, it 
suffices to show that the condition (2.22) holds on this slice, where A(t) is the large N instanton 
action corresponding to the first non-trivial instanton sector with N 2 = 1. This action is given 
in (3.40), and it resums quantum fluctuations around the classical instanton solution. At small 
t it should have the general structure [23] 
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where the prefactor is the one-loop fluctuation around the classical instanton with action c\ = 
ir 2 /2, c%N is the number of zero modes at large N, and C2 is a numerical constant. Therefore, 
the expansion of A(t) around t = should be of the form 

A(t) = c 1 -c 3 t\og(^)+0(t). (4.8) 



t 

Indeed, we find from the exact expression (3.40) that 

A(t) = ^-tlog^j -t + 0(t 3 ). (4.9) 
This large N instanton action is real for t G [0, 00) and it decreases monotonically from 

7T 2 

A(0) = — to lim A(t) = 0. (4.10) 

2 t— >oo 

In particular it does not vanish for any finite t. This can be seen in Fig. 2, where the continuous 
line represents A(t)/t as a function of t/ir. As we explain in the Appendix, this is equivalent to 
the absence of phase transition for q-deformed 2d Yang-Mills with p = 2 [12, 5, 28]. We conclude 
that, for real g s , the large N asymptotics of the partition function of Chern-Simons theory on 
L(2, 1) is simply given by the 1/N expansion of Chern-Simons theory on § 3 , 

00 

\ogZ L{2tl) {N,g s ) = Y,gl 9 ~ 2 Ff{t). (4.11) 

3=0 

In other words, the relevant saddle geometry along this direction in moduli space is just the 
resolved conifold geometry. 

Of course, the result (4.11) has an infinite number of exponentially suppressed instanton 
corrections of the form (2.11). From the point of view of classical asymptotic analysis, these 
subleading saddles are not taken into account, but of course they are important thanks, among 
other things, to the Stokes phenomenon that we will now uncover. 




Figure 1: In this figure we plot log \Z(N, N2,g s )\ as a function of iVjj = 0, • • • , N, for N = 12 and for 
g s = ±7r/24 respectively on the left and on the right. The dominant configurations are (N*, JV^ ) = (N, 0) 
and (0,N), respectively, in agreement with the symmetry (4.3). 

The above analysis can be confirmed numerically for small values of N. In general, a simple 
way to estimate the leading saddle (N* , iV| ) for finite is to calculate the partition functions 
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Figure 2: In this figure, the crosses represent the right hand side of (4.12), for N — 24 and different 
values of t, while the continuous line represents the value of A(t)/t obtained from (3.40). 



Z(N, N2, g s ) for fixed values of N 2 , and see which one is the largest in absolute value, and 
therefore gives the dominant contribution in (4.2). In Fig. 1 we plot log \Z(N, N2,g s )\, for 9s real 
and positive. As expected, the dominant configuration is (iVf, iV|) = (N,0) ((Ni[,N%) = (0,N) 
for g s real and negative). 

We can also confirm numerically the value of the instanton action (3.40). From (2.11) we 
find the large N asymptotics 

1 f Z(N,N 2 = l,g 8 ) \ A(t) (1\ 

In Fig. 2 we plot the values of the r.h.s. for N = 24 and different values of t = Ng s , and 
we compare them to the value of A(t)/t given by (3.40). As we can see, the agreement is 
excellent, and it can be further improved by extracting the subdominant tails in (4.12) through 
the technique of Richardson transforms. 

4.2 Imaginary g s : the interior expansion 

Let us analyze now the large N asymptotics along the direction Arg(g s ) = tt/2. The boundary 
point (iV, 0) can not be a solution any longer, since (2.22) is no longer verified (at least for small 
t). This is easy to understand: for imaginary g s , (3.11) is purely imaginary and the non-trivial 
instanton sectors are not suppressed anymore. We then have to look for saddles satisfying the 
condition (2.24), which in this case reads: 



Rc 



g s V dti dti 



0. (4.13) 
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It is easy to see that the real part of i*b(*l ? ^2) is symmetric under the exchange of t\ and t 2 , 
therefore the configuration 

*i=*2 = |, (4.14) 

which is an interior point, is a saddle. We claim that this gives in fact the dominant saddle 
when g s is imaginary. Of course, there could be other saddles which actually dominate the 
asymptotics. A first indication that this is not the case, and that (4.14) is the relevant saddle, 
comes from the numerical analysis of Z(N, N 2 , g s ) for small N. As we see in Fig. 3, for N even 
the largest contribution comes from N 2 = N/2, while for N odd the dominant contributions have 
N 2 = (N±l)/2. 



Figure 3: In this figure we plot log | Z(N, N2,g s )\ as a function of N2 = 0, • • • , N, on the left for N = 12 
and g s — 7ri/48, on the right for N — 13 and g s = ni/52. The largest value is obtained when N2 = N/2 
for N even and when N = (N ± l)/2 for N odd. 

According to the general discussion in section 2, since the saddle (4.14) is an interior point, 
the large N asymptotics of the partition function should be given by (2.17). In our case, we have 

1 u tt 2 irik 

6 = -, » = v = 0, - = — = - (4.15) 

2 g s Zg s 2 

The resulting theta function depends on the value of N modulo two. If N is even, Ne is an 
integer which can be reabsorbed in a shift of n £ Z, which is summed over. Therefore, (2.18) is 
given by 

(k7T \ 
■Kirn 2 + i — n J . (4-16) 
/ 

Notice that, since g s is imaginary, k is real, and the second term leads to an oscillatory behavior. 
If N is odd, Ne is a half-integer. We can absorb its integer part by a shift of n and we obtain for 
(2.18) 

v / k-K \ 

■& = 22 exp 7rirn2 + i ^ _n • ( 4 - 17 ) 



We then claim that, when g s is imaginary, the asymptotic behavior of F(N, g s ) = log Z(N, g s ), 
the free energy of Chern-Simons theory on L(2,l), at large N and fixed 't Hooft coupling, is 
given by the logarithm of (2.17). At leading order we have 

F(N, g s ) = g 2 s F f |, |) + F 1 ^ + logG ,o(n, r) + • • • (4.18) 
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where t = Ng s . Notice that 

1. The perturbative free energies F g (t/2,t/2) are the genus g free energies of topological string 
theory on Fo on the slice t\ = t<i = t/2. For genus zero, it is given by the explicit expression 
(3.42). The higher genus amplitudes are essentially quasi- modular forms of the modular 
parameter r [2, 25] which is given in (3.43) by an explicit function of 

2vriiV 

n = k + N' (4 ' 19) 

where we used the original Chern-Simons parameters. In particular, F\ is given in (3.35). 

2. The asymptotics depends on the parity of N through the theta function, which is given by 
(4.16), (4.17) for N even and odd, respectively. 

In the quantum Chern-Simons gauge theory, k = k + N is an integer, and g s is purely 
imaginary. Therefore, the above asymptotics is in principle the relevant one to understand the 
large N 't Hooft limit of the quantum Chern-Simons invariant of L(2, 1). 
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Figure 4: On the left column we plot the sequence Hef e (£,z) as well as its Richardson transforms 
Re/i fe ^, z) for k = 1,2, 3, and for £ = 2- • • 14. On the right column we plot the sequence Re/ (£, z) as 
well as its Richardson transforms Re fo h \i, z) for k = 1, 2, 3 for £ = 2 • • • 14. The plots on the top are for 
z = 1/4, and the plots of the bottom are for z — 1/8. The dashed lines in the plots show the expected 
value -ReF (z)/(7rz) 2 . 

Our claim about the asymptotics follows from the identification of (4.14) as the relevant 
saddle point, and from the general theory reviewed in section 2 concerning instanton corrections. 
We now provide detailed numerical evidence for our claim. Instead of working with N and g s we 
will use N and z, which is defined by 

N 

z = T . (4.20) 
k 
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Let us now define the sequences of complex numbers obtained from the non-perturbative partition 
function with N = 2£, 21 + 1 respectively: 

m Z) = j^f l0g Z(2e > Z) > m Z) = (21+1)2 log Z{21 + *)" ( 421 ) 
According to our claim, these sequences have the large £ asymptotics 



f e (£, Z ) = + (f 1 (z) + logMz)) + o(r 4 ), 



(4.22) 



where Fq : \{z) are the genus 0, 1 free energies evaluated on the slice (4.14). Notice that the theta 
functions have in general oscillatory behavior as a function of £, but this can be fixed by choosing 
special values of k. For example, if k = mod 4 we have: 

e (z) = 3 (r), #o{z) = ±Mt), (4.23) 

where the sign depends on the value of k mod 8. In this way, the resulting large i asymptotics 
is simply a series in inverse powers of £, and we can use standard techniques of accelerated 
convergence to compare the actual values of / eo (^; z )i computed from the exact partition function 
of Chern-Simons theory for low £, to the predicted asymptotics. 




Figure 5: The sequence ReO eo (z) for z = 1/4, 1/8, together with its Richardson transforms ReO^o (z) 
for k = 1, 2, 3. The dashed lines in both plots show the expected value Re(log , )? e (z) — \og'd (z)). 



In practice, the computation of (3.5) involves generating a large number of configurations for 
the vectors n, and we obtained their numerical values for N = 1, • • • ,30 and for different values 
of z. This means that we obtained the values of the sequences f e ,o(£, z) up to £ = 14. Once these 
sequences are computed, we obtain their Richardson transforms (see for example [6]) 



E 

fc>0 



f e , (£ + k,z)(£ + kY(-l) k+r 
k\{r-k)\ 



(4.24) 



in order to accelerate the convergence. In Fig. 4 we plot the sequences Re f e:0 (£, z) as well as 

their Richardson transforms Re/i^j }(£, z) for k = 1,2,3 for different values of z together with 
their expected value — Re Fq(z)/(ttz) 2 . As we can see, the agreement is excellent. 
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Figure 6: On the left column we plot the sequence Re/ ei i(£, z) as well as its Richardson transforms 
Re/g^ i(£, z) for k = 1,2, 3, and for £ = 2- ■ ■ 14. On the right column we plot the sequence Re/ 0j i(^, z) as 

well as its Richardson transforms Re f^ k }(£, z), again for k = 1,2,3 and for I = 2 • • • 14. The plots on the 
top are for z = 1/4, and the plots of the bottom are for z — 1/8. The dashed line shows the expected 
value Re(-Fi(2:) + \ogd e , (z)). 



In order to test the subleading behavior, we define the following sequences: 

Q eo (£,z) = (2£) 2 (f e (£,z)-f (£,z)), 



f eA (£,z) = (2£) 2 [f e (£,z) + 



[TTZ 



(4.25) 



f oA (£,z) = (2£ + l) 2 (fo(£,z) + 



Fo(z) 
(ttz) 2 



with the expected asymptotic behavior: 



lim & eo (£,z) = logi9 e (z) - logtf (z) 

l— >oo 

lim f e>1 (£, z) = F 1 (z)+ log $ e (z), 
lim f 0tl {£,z) = F 1 (z) + log# (z). 



(4.26) 



As we can see in Fig. 6, the agreement between the analytic prediction and the actual asymptotic 
behavior is again extremely good. This confirms our identification of the saddle point (4.14), 
as well as the formalism of [10, 17, 18] to incorporate instanton corrections in the oscillatory 
case. Although these corrections are subleading in the case of the free energy, they can appear 
at leading order when studying expectation values of operators like Wilson loops, as noticed in 
[10] in the context of matrix models. 
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4.3 The phase diagram for complex g s 



A 

t\ 




o e tt/2 

Figure 7: The phase diagram describing the large N limit of Chern-Simons theory on L(2, 1) for < 8 < 
7r/2, where 8 — Arg(<7 s ). There are two phases distinguished by the equilibrium value of t 2 and separated 
by the critical curve (4.30). In the region to, — 0, the 1/iV asympotics is given by a genus expansion, 
and corrections due to neighbouring geometries are exponentially suppressed. In the region where t 2 ^ 0, 
neighboring geometries correct the genus expansion with oscillatory terms. Notice that, in this region, the 
equilibrium value of t 2 changes as we change t. The line 8 = tt/2, where t 2 — t/2, can be regarded as an 
anti-Stokes line. 

The qualitative features of the phase diagram for < Arg(g s ) < tt/2 can be understood by 
looking at the function (3.40). Let us denote 

9 = Avg(g s ) (4.27) 

so that 

t = \t\e ie . (4.28) 

If Re(A(t) / g s ) > 0, when t, g s are now complex, the dominant saddle will be t\ = t, t% = 0, as 
for 9 = 0. For 7r/2 > 9 > 0, the function 

f (\t\)=Re(e- i9 A(e w \t\)) (4.29) 

vanishes for a certain critical value 

\t\ = \t\ c (9), (4.30) 

which can be easily determined numerically. The curve (4.30) separates two different regions in 
the complex t plane, as shown in Fig. 7. In the region under this curve, where 

< \t\ < |i| c (0), (4.31) 

the dominant saddle is t\ = t, ti = 0. In the region above this curve, we expect the dominant 
saddle to have t2 / 0. 
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Figure 8: Plot of (4.32) (blue), calculated for N = 22, and the real part of the instanton action (4.29) 
(red), as a function of \t\/ir, for 6 = 97r/20 , 77r/20 and tt/4, from left to right. The function (4.29) 
decreases monotonically until it becomes zero at |t| c (0)- At this point, the average (4.32) gets triggered 
and it increases with \t\. 



In order to verify this, as well as the location of the phase boundary, we have computed 
numerically the average filling fraction at finite N 

for fixed 0, as a function of \t\. This gives a finite N approximation to the value of \t2/t\ at 
the dominant saddle. In Fig. 8 we plot the value of (4.32) and the real part of the instanton 
action (4.29), for three different values of 0, and as a function of \t\/ir. Within the limits of 
the numerical approximation, we clearly see that N2/N starts developing an expectation value 
precisely when we reach (4.30). As shown in Fig. 9, the limiting value of \t2/t\ as \t\ increases 
seems to be 1/2. 




Figure 9: Plot of (4.32) as a function of \t\/n, for N = 22 and = 7tt/20. 



- 23 - 



arg(# s ) = 



Z{N,g s ) 




t 





+ 



+ ••• 



arg(^ s ) = tt/2 




Z{N,g s ) 



+ ••• 



Figure 10: The Stokes phenomenon for the large N duality of L(2, 1). When g s is real and positive, 
the saddle geometry giving the dominating contribution to the non-perturbative perturbative function is 
simply the resolved conifold. The corrections, which are exponentially suppressed, can be interpreted as 
background geometries in which t 2 = lg s is quantized in units of g s . When g s is purely imaginary, the 
saddle geometry is a local P 1 x P 1 with t\ — t 2 . Corrections to the genus expansion contribute to the 
large N asymptotics. They can be regarded as a sum of infinite geometries where i 2 = n g s is quantized, 
and their contribution adds up to a theta function. In terms of the large radius Kahler parameters, these 
geometries have fixed T\ — T 2 = t. 

4.4 Stokes phenomenon and target geometries 

We can now summarize the most important qualitative results of the analysis of the phase 
diagram. 

First of all, for each complex value of t, the large N asymptotics of the non-perturbative 
free energy is dominated by a fixed filling fraction or background (t\,t\). Such a background 
corresponds to a fixed target geometry, i.e. to a point in the Calabi-Yau moduli space. However, 
the dominant background changes as we vary t. 

The second important point is that the structure of the large N asymptotics depend on the 
nature of the dominant saddle point. If the dominant saddle point is on the boundary, i.e. one 
of the filling fractions vanishes, the 1/N asymptotics is given by a conventional genus expansion. 
There are non-perturbative effects due to the contributions of neighbouring geometries, as shown 
in (2.11), but they are exponentially suppressed. However, if the dominant saddle point is an 
interior point, the asymptotics has an oscillatory behavior and it involves theta functions, as 
detailed in (2.17). In particular, the large N asymptotics is no longer given by a genus expansion 
around a fixed geometry, and corrections due to neighboring geometries are crucial already at 
the next-to-leading order. We hope that our detailed numerical analysis in the case 6 = tt/2 
has convinced the reader that these effects are crucial in order to obtain the correct large N 
asymptotics of the gauge theory partition function. 

It is interesting to notice that the non-perturbative corrections are due to fluctuating target 
geometries whose size is quantized in units of g s , see Fig. 10. In the case of the boundary 
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expansion, an £-instanton correction corresponds to a geometry with ti = lg s , £ = 1, 2, • • • . In 
the interior case, the theta function is a discrete sum over an infinite number of geometries with 
h = (n — N/2)g s , n G Z. Notice that the saddle geometries for arg(<7 s ) = tt/2 correspond to the 
region in moduli space where z±,Z2 — > oo with Z1/Z2 fixed and finite, therefore they are not in 
the large radius phase z\, zi — > 0. 

Finally, we point out that the change of saddle geometry as we change the complex 't Hooft 
parameter can be regarded as a generalized Stokes phenomenon. In fact, in the limit of vanishing 't 
Hooft coupling, the asymptotics changes discontinuously as we change the argument of the string 
coupling constant: we have a sudden jump from a dominant saddle at t2 = to a dominant saddle 
at t2/t = 1/2. The jump in the asymptotics has the same origin as in the Stokes phenomenon: 
saddles which were subleading due to Re(A/g s ) > are no longer suppressed exponentially, and 
they lead to an oscillatory asymptotics in the 8 = tt/2 direction. In fact, the line 6 = tt/2 plays 
the role of an anti-Stokes line. These analogies can be made more precise by studying a model 
where the connection to the standard Stokes phenomenon is manifest and can be followed in 
detail, namely the cubic matrix model. We will show in the next section that the phase diagram 
for this model is very similar to the one we have just found, and we will able to confirm all the 
analogies that we have stated. 

It is important to notice that the connection to the standard Stokes phenomenon is only 
strictly true when t — > 0. For finite \t\ the Stokes phenomenon is smoothed out, in the sense that 
for fixed \t\, as we increase 6 from 9 = 0, the asymptotics changes smoothly from a phase with 
t2 = to a phase with t2 7^ 0. The phase boundary (4.30) corresponds then to a second-order 
phase transition with t2 as its order parameter. 

4.5 Background independence 

Our results on the phase structure of this model shed some light on the issue of background 
independence in topological string theory. As it is well-known, the genus g free energies in 
topological string theory, F g (ti), are quasi-modular forms under the action of the symplectic 
group [2]. This can be regarded as a consequence of the holomorphic anomaly of [8], and it 
means in particular that the F g (ti) are not well-defined functions on the moduli space. 

In the models we are studying, thanks to the large N duality, there is a quantity which is well- 
defined non-perturbatively, namely the Chern-Simons theory partition function, which depends 
on g s and the rank N of the gauge group. Notice that, as needed for consistency, the total 't Hooft 
parameter t is a symplectic invariant in the dual topological string. An important question is how 
this non-perturbative quantity is related to the perturbative, background dependent topological 
string free energies. In this paper we have provided a detailed answer to this question in the 
particular example of L(2, 1) and its string dual. 

First of all, at finite N the non-perturbative partition function (2.7) is a sum over all possible 
backgrounds N1+N2 = N. In the gauge theory description, these backgrounds correspond to sad- 
dle points or flat connections. In the matrix model description, they correspond to different filling 
fractions. In the language of statistical mechanics, we can say that the background-independent, 
non-perturbative partition function has to be calculated in a grand-canonical ensemble where 
N2 is a fluctuating variable. This is manifest in (2.7), where the parameters Q play the role of 
fugacities. On the other hand, the background-dependent topological string partition function is 
computed in a canonical ensemble where N2 is fixed. 

In the thermodynamic limit (i.e. at large N), as it is well-known in statistical mechanics, 
the sum over backgrounds is strongly peaked around a particular one, and both ensembles are 
equivalent. This particular background is the dominant saddle which we have determined in 
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various situations. It depends on the gauge theory parameters ./V and g s , and it changes as we 
move on the complex t plane, displaying the rich phase structure that we have analyzed. At finite 
N, however, the results of the two ensembles are different, and they differ precisely in the terms 
which go beyond the genus expansion, i.e. in the contributions of large iV multi-instantons. 

The resulting non-perturbative picture is very different from the picture that one obtains 
in perturbative topological string theory: the background-dependent quantities are emergent 
quantities and they only make sense in the large N limit, since they are defined through an 
asymptotic expansion. 

5. Stokes phenomenon and matrix models 




Figure 11: The integration contour C for (5.2), which leads to a convergent integral. 

In order to have a better understanding of the nature of the line Avg(g s ) = tt/2 in the phase 
diagram of Chern-Simons theory, we will make a detailed analysis of the Stokes phenomenon in a 
simple matrix model. We will focus on the quintessential example of Stokes phenomenon, namely 
the Airy function, and we will study its fate as we promote the Airy integral to a matrix integral. 
Of course, the study of matrix models in the double-scaling limit produces differential equations 
which display the so-called non-linear Stokes phenomenon [39, 21]. Here we are rather interested 
in recovering the standard, linear Stokes phenomenon in the context of matrix integrals. 

5.1 Review of the Stokes phenomenon for the Airy function 

Let us consider the potential 

3 

V(x) = -e iK x + \ (5.1) 
where k £ [0, 2n] . The one-dimensional integral 
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Figure 12: Saddle-point analysis of the asymptotics of the function (5.2), closely related to the Airy 
function Ai(£). Full lines (in red) are Stokes lines, while dashed lines (in blue) are anti-Stokes lines. 
On the Stokes lines K — ±2tt/3, a second saddle appears in the integration contour. This saddle is 
subdominant when 2tt/3 < \n\ < n and does not contribute to classical asymptotics. However, at k = n, 
the saddle is not subdominant anymore and leads to an oscillatory asymptotics. 



where C is the path shown in Fig. 11, and g s is real and positive, is essentially the Airy function 
Ai(w), where 

w = 5 7 2/3 e iK . (5.3) 

The small g s asymptotics of this integral is given by the large w asymptotics of the Airy function. 

The standard analysis of this asymptotics is as follows (see [38] for a very nice discussion). There 
are two saddle points 

X L,R = tC !/ 2 (54) 

where we have introduced the variable 

C = e iK . (5.5) 

These saddles have actions 

--V{x^) = ^e'\ (5.6) 

For 

- T <«< T (5.7) 

the path C can be deformed into a path of steepest descent through the saddle point at x L . 
When k = 27r/3, the steepest descent path coming from the saddle at x L runs right into the 
other saddle point. At this angle we have 

lm{V(x L )) = Im(T/(x R )) (5.8) 

and the corresponding direction is called a Stokes line. This is the place where the second saddle 
x R might start contributing to the integral. In fact, for 

2tt . , 

— < \k\ < vr (5.9) 

O 
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the contour C gets deformed into a steepest descent path passing through x L together with a 
steepest descent path passing through x R . However, in this range the latter gives an exponentially 
suppressed contribution, and the classical asymptotics is just given by the contribution from x L . 
This leads to 

Zl ~e- 2w ' i/2 / 3 (5.10) 

where w is the variable introduced in (5.3). When k = tt the asymptotics is different, since both 
saddles have the same real part 

Re(V{x L )) = Re(V(x K )). (5.11) 

A line where this occurs is called an anti-Stokes line. Along this line both saddles contribute to 
the asymptotics, and we have 

Zi ~cos(jM 3/2 - J). (5.12) 

The fact that different asymptotic formulae hold on different directions for the same analytic 
function is called the Stokes phenomenon. From the point of saddle-point analysis, what is 
happening is that the saddle point which appeared on the Stokes lines, at k = ±2-7r/3, is no 
longer subdominant at k = tt, and it has to be included in the asymptotics. The saddle-point 
analysis is summarized in Fig. 12. We should mention that it is important sometimes to take 
into account the subleading exponentials, even in the region k / tt. This is very clearly shown 
in the analysis of Berry in [7], and it has led to the development of non-classical asymptotic 
theories. The most sophisticated of these, resurgent analysis, is applied to the Airy function in 
for example [15]. 

5.2 The Stokes phenomenon and the cubic matrix model 

Let us now promote (5.2) to a full matrix integral, 

„ N , 

Z C np (N, 9s ,n) = P^^) 6 ^^ 1 ^ 1 ' ( 5 - 13 ) 

where all the eigenvalues are integrated along the contour C. We can now study its small g s 
asymptotics for t = g s N fixed. Of course, this is nothing but the 't Hooft expansion of this 
matrix model. The general saddle point is a two-cut configuration, labelled by (N\,N2), where 
N±, N2 are the number of eigenvalues near the critical points x L ' R , respectively. Before starting 
the analysis, let us notice that for small t the saddle point structure should be the same as for 
the Airy function. The reason is that, for small t, the Vandermonde repulsion among eigenvalues 
is suppressed w.r.t. the potential, and the model becomes just N copies of the one-dimensional 
integral. 

Therefore, at least for small t, we expect the following phase structure. For \n\ < 2tt/3, the 
dominant saddle is a one-cut configuration where all the eigenvalues sit near x = -C 1/2 . This 
is the boundary saddle point (N,0). It is a one-cut configuration which can be analyzed with 
standard techniques. The endpoints of the cut (a, b) are determined by the equations 

x (x 2 -Q = t, 6 2 = 2(C - xl) (5.14) 

where 

a = —xq + 5, b = —xq — 5 (5.15) 
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and we choose the root xq = C}^ 2 + ■ ■ ■ . The spectral curve is 



y (x) = (x- x )^Jx 2 + 2xx + 3x 2 - 2( (5.16) 
and the effective potential reads 



1 



Vce(x) = - (x(x - s ) - 2C) V x2 + 2xx o + 3x 2 - 2C 

- 2x (x 2 - C) log ^x + x + y / x 2 + 2xx + 3xg-2C^ . 



(5.17) 



For 27r/3 < \k\ < ir, the second saddle should start contributing, but at least at small t it should 
be exponentially suppressed. In other words, the saddle configuration for the matrix integral in 
this region should still be the boundary saddle (N,0), but there will be corrections of the form 
(2.11). The instanton action is given by 

A = V cS (x ) - V cS (b). (5.18) 

The small t expansion of this action is 

A = ~- 3 C 3/2 -t + tlogt- 3Hog(2C 1 / 2 ) - ^L_t 2 + ^ + ■■■ (5.19) 

Notice that, as t — > 0, we find 

x -»• x R , (5.20) 

which is the other saddle, and 

A^-^ 2 (5-21) 

which is the difference between the actions of the two saddles (5.6). This is in agreement with 
the expectation that for small t the saddle-point structure of the matrix model is just the one 
coming from the Airy function. The real part of the instanton action vanishes at a critical value 
of t c (/c), which depends on the value of k. For example, for k = 27r/3, we have 

'°=iM- (5 - 22) 

This is precisely the critical point leading to pure 2d gravity. For t > t c (/c), a phase transition 
occurs to a new phase, in which we expect generically N% ^ 0. 

When k = 7r, i.e. on the anti-Stokes line, we expect that the standard Stokes phenomenon 
controls the behaviour of the matrix integral for small t. Indeed, this is the case. The real part 
of the action A(t) vanishes at t = 0, and it is actually negative for small t. We then have to look 
for a new saddle point by analyzing the behavior of Fq. The genus zero free energy of the cubic 
matrix model has been analyzed in many places [11, 30]. It has the structure 

Fo(h,t 2 ) = F^{h,t 2 ) + F^(h,t 2 ), (5.23) 

where 

Fo P (h,t 2 ) = it? (in ti-|)+ ^(ln t 2 - 3 -)- \{t\ + t 2 ) log(2 C 1 / 2 ), (5.24) 
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Figure 13: The phase diagram of the cubic matrix model (5.13) as a function of 27r/3 < k < it and 
t = g s N (which is taken to be real), for small t. The Stokes and anti-Stokes lines of the N = 1 case, which 
occur for k — 2n/3 and k — 7r, respectively, determine to a large extent the phase diagram for small t. On 
the anti-Stokes line, the saddle value of t 2 is t/2, and the classical 1/N asymptotics requires corrections 
beyond the genus expansion. 



and 



t 2 ) = -^C 3/2 (h - t 2 ) + 2ht 2 log(2C 1/2 ) + (--t? - 5t|ti + 5t 2 t? + -i| 



1 /8 4 91 o » a 91 o 8 „ 



(5.25) 



On the anti-Stokes line arg(«) = tt, C 1//2 i s purely imaginary, and it is easy to see from the 
structure of the genus zero free energy that its real part is symmetric in t±, t 2 - Therefore, 

ti = (5.26) 

solves the saddle-point equation 

Re ^ = (5.27) 

OS 

at least for small t (more precisely, (5.26) provides a solution in the domain of analiticity of the 
convergent expansion (5.25)). Therefore, on the anti-Stokes line and for t small enough, the 
partition function of the cubic matrix model (5.13) is given by an expansion of the form (2.17). 
This type of expansion is then the natural generalization to matrix models of the oscillatory 
behavior along an anti-Stokes line. 

The phase diagram of the cubic matrix model, again at small t, is represented in Fig. 13. 
We see that it is very similar to the phase diagram of the large N Chern-Simons theory on a 
lens space that we analyzed in the previous section. This confirms that the line 9 = vr/2, where 
the vacuum occurs at t\ = t 2 = t/2, plays the role of an anti-Stokes line. Strictly speaking, 
however, the discontinuous jump in the asymptotics typical of the Stokes phenomenon occurs 
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only at t = 0. For t 7^ 0, as we change the angle we encounter the region where ti goes from a zero 
value to a nonzero value in a smooth way. We conclude that turning on the 't Hooft parameter 
smooths out the Stokes discontinuity, and we pass from a first-order phase transition at t = to 
a second-order phase transition for t > 0, where ti plays the role of an order parameter. 

As t increases along the different directions, further phase transitions can occur which we 
have not explored. For example, along the k = tt direction, it follows from the results of [9] that 
there is a transition at t > in which the two cuts merge into a trivalent graph. Interestingly, 
such transitions along the anti-Stokes line are absent in the lens space partition function. 

6. Conclusions 

In this paper we have analyzed non-perturbative aspects of the large N expansion in a theory 
with a rich structure of instanton sectors and with a large N dual. There are various important 
conclusions that emerge from our analysis. First of all, we have seen that, for each choice of 
gauge theory parameters g s , N, the large N asymptotics is dominated by a single instanton 
sector. This sector corresponds to a fixed geometry in the moduli space of the dual Calabi-Yau 
geometry, which can be characterized by the Boutroux condition (2.25). Second, there are large 
N phase transitions as we move in the complex t plane, and different instanton sectors dominate 
the asymptotics for different values of t. These transitions can be regarded as deformations of 
the standard Stokes phenomenon. Third, the correct large N asymptotics goes beyond the genus 
expansion, and subleading saddles must be incorporated. 
We conclude with a list of issues left open by our work: 

1. The genus g topological string free energies appear in the large N asymptotics of the 
"canonical" partition function Z(N, N2,g s ) with N2 fixed. They are known to transform as 
quasi- modular forms when moving from one patch of the moduli space to another [2]. At 
the level of the gauge theory, these transformations should correspond to phase transitions 
within the same instanton sector, which we have not studied in detail in this paper, since 
we have focused on phase transitions among different instanton sectors. On the other hand, 
the matrix integral (3.8) gives a non-perturbative definition of this "canonical" partition 
function which is valid for any and any complex g s . It would be important to understand 
the precise connection between the modular properties of the topological string amplitudes 
and the large TV" asympotics of the quantity (3.8), which is background independent and 
clearly well-defined in all the moduli space. The results obtained in this paper seem to indi- 
cate that modularity can be understood in terms of phase transitions or generalized Stokes 
phenomena, and that background dependence in these models is ultimately an artifact of 
the asymptotic large N expansion. 

2. Non-perturbative effects of the form e~ A ^ 9s in gauge theories and matrix models have been 
identified in terms of D-branes in their string duals. In particular, the source of these effects 
has been shown to be Liouville branes in minimal strings, and more recently toric branes 
in topological strings [33, 35, 43]. It seems clear that the subleading saddles appearing 
in the oscillatory asymptotics (2.17) can be also interpreted in terms of D-branes, and 
in the example of local Fo these arc also toric D-branes in this geometry. It would be 
very interesting to make more precise the D-brane interpretation of the non-perturbative 
corrections appearing in (2.17). 

3. Even though in this paper we focused on a simple toy model of gauge/string duality, we 
believe that some of the effects we studied in this work are likely to play a role in other 
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large N dualities. The underlying reason for the phenomena explored in this paper is the 
existence of nontrivial instanton sectors which compete among them at large N. In some 
cases we have a single sector dominating the large TV asymptotics, but in other cases, along 
generalized anti-Stokes lines, the standard large ./V asymptotics needs to be corrected, as 
first pointed out in [10]. This scenario seems to be quite general and it is likely to appear 
in other models, in particular in AdS/CFT dualities. For example, the recent work [29] 
indicates that matrix integrals closely related to (3.9) are relevant in the study of large N 
dualities for ABJM theories. 
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A. Some aspects of topological string theory on the local P 1 x P 1 geometry 

In this Appendix we give some further information on topological string theory on the local 
P 1 x P 1 geometry, and in particular we justify some of the statements made in the bulk of the 
text. 

The Yukawa couplings of this Calabi-Yau, when written in terms of the "bare" coordinates 
z±,Z2 appearing in (3.20), are [3, 25] 

_ (l-4z 2 ) 2 -16zi(l + zi) 
t^m — 



Cli2 — 



4zfA 
16zf - (1 - Az 2 f 



4zfz 2 A 

(A.l) 

164 -(l-4.i) 2 
° 122 ~ 4^|A ' 

(1_ 4zi) 2 -16z 2 (l + * 2 ) 

^222 — 



4zfA 
where 

A = 1 - 8(zi + z 2 ) + 16(zi - z 2 f. (A.2) 

Using the Yukawa couplings (A.l), it is possible to find a closed form for the prepotential 
when x 2 = 0. The triple derivative w.r.t x\ is given by 



2 2 1 — X\ 

= (CiuJ 1 + 3Cii2 J\ J 2 + "2>C\ 22 3\J 2 + C 222 J 2 ) „ = - — — o, (A. 3) 

x 2 =o v /x 2 -o 4xi (xi - 1) 



where 

dz\ x\ — 2 dz 2 
dx\ x\x\ ' 2 dx\ x\x\ 



J i - ~ .-3^.2 ' J 2 - j^zr - ( A - 4 ) 



Since 



Sxi o"i = ~, , (A.5) 

1 — X\ 
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we obtain a closed form for the triple derivative w.r.t. the flat coordinate a±: 

1 + e" ai 



a 



X2 



=0 4(l-e- CT i)' 



which can be integrated to give: 



^o(^i,0) = --Li 3 (e- <Tl ) + 



(A.6) 



(A.7) 



The dots indicate a degree three polynomial. This polynomial can be fixed by comparing to the 
perturbative results, and one finally obtains (3.42). 

It is also possible to derive closed expressions for d a2 T a2 on the slice x 2 = 0. To do this, one 
first writes the Picard-Fuchs system (3.21) in terms of the coordinates xi i2 , 



d = - (8 - 8xi + x\) x 2 8 X2 + t (-4 + (-2 + xi) 2 xfj d 2 



x 2 



4 — i/ 4 
+ (-1 + xi) x\d xi - xi (2 - 3xi + x\) x 2 d Xl d X2 + (-1 + xi) 2 x\d 2 



(A.8) 



£ 2 = (2 - xi) x 2 d X2 + (-1 + (1 - xi) x\) dl 2 - x\d xi 
+ 2 (-1 + xi) xix 2 d Xl d X2 + (1 - xi) x\d 2 Xl , 
and looks for a solution £1,2/ = in the form 

/(xi,x 2 ) = x 2 G(xi)+C(x|). 
The PF system leads to the following equation for G(x): 

-Ax -A) G{x) + x G'{x) - (1 - x)x 2 G"(x) = 0, 



(A.9) 



(A.10) 



whose general solution is a linear combination of X-fsT(x) and xK(l — x). By comparing to the 
first terms in (3.26) and using (3.30), one immediately obtains 



therefore 



da 2 
dx 2 

d 2 Fn 



= -xi-fT(xi), 

x 2 =0 71 



dx 2 



ds 



= -x\K(l - xi), 



X2=0 



ds 2 



= 4 



x 2 =0 



da 2 

dx 2 



-1 



3:2=0, 



dx 2 V ds ) X2=0 n K( Xl ) 



(A.11) 



(A.12) 



and (3.43) follows. Equivalently, one can derive (A. 11, A. 12) by solving the recursion (3.27) for 
Cm,i and dm^i, 

r(m-l/2) 2 

-, m > 1, 

(A.13) 



Cm,l 



^r(m) 2 ' - ' 

dm,l = -2(Y> m - ^m-l/2) c m,l, 



and using the series expansions 



K(x) 



7T 



00 , 

El 



k=0 



2 fc 
X K 



- x ) = -- \og{x)K{x) + 2^ ^2 J fe (fe!)2 ~ ^+V2)- 



(A.14) 



- 33 - 



Finally, we calculate the instanton action (3.40) in the slice t 2 = 0. This can be in principle 
derived from the above geometry, but in fact we can use a computation already done in [5] . As in 
[35, 36], the one-instanton free energy is given by the ratio of the one instanton partition function 
Z*jp = Z(N,N2 = 1,5s) and the zero instanton partition function zff = Z(N,N2 = 0,g s ), 



7(1) 

F (i) = £jv_ 
7 (o) 

N 



~ e 



-A{t)/g s 



(A.15) 



We can evaluate zff by taking (Aq, Af 2 ) = (N — 1, 1) in the matrix integral (3.9), as in [35], and 
we find 



7(1) 



N ~ 1 1 2vr 




2 sinh 



S + 7T1 — Xj 



(A.16) 



i=2 



N-l 



where pN is the Weyl vector in the U{N) theory, and the vev is calculated in the matrix model 
defined by zff. At small g s , the integral over s can be calculated at the saddle s = 0, and the 
v.e.v. can be calculated in the planar limit. Taking into account that 



7(0) 

7 (0) 

N 



cxp 



.9. 



12 



(A.17) 



we find, 



where 



^) = y + S t (if(t)-^ 



;ydup(u)log(^4cosh 2 ^, (A.18) 



p(u) 



TTt 



tan 



-1 



cosh 2 (f) 



cosh (^ 



(A.19) 



is the density of eigenvalues of the S 3 Chern-Simons matrix model. 

Fortunately we do not have to evaluate the integral in (A.18), since the result is already 
available in the literature. A closely related integral has been studied in [5, 12, 28] in the context 
of the g-deformed two-dimensional Yang-Mills theory on the two-sphere, and it calculates the 
large N instanton action 7(^4, p). Its explicit expression can be found in equation (4.22) of [5]. 
In particular, for p = 2 and A = — At, (A.18) equals 7 (.A, p)/4, and using the result of [5] one 
finds 



Ait) 



7T 



2Li 2 (l - e"* /2 ) + 2Li 2 (e"* /2 ) + t log(l - e" 



*/2\ 



Li 2 (e" 



= 2(Li 2 (e-'/ 2 )-Li 2 (-e-*/ 2 )). 



(A.20) 



For p = 2, the g-deformed Yang-Mills theory does not undergo a large N phase transition, as 
shown in [5, 12, 28]. In our context this means that there are no large phase transitions when 
t is real. 
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